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ABSTRACT 

Physical data intended to represent the variation of a function of a single 
variable may be actually mean values of the function over intervals of the argu- 
ment. Formulas for approximating to the values of the function from such 
data and examples of their use are presented in this paper. These formulas have 
been applied to problems of strain distribution in the Engineering Mechanics 
Section of the National Bureau of Standards. 
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I. INTRODUCTION 

Methods of measurement to determine the variation of a function 
with its argument are frequently such that only the mean values 
of the function over intervals of the argument are obtained. For 
example, a strain gage measures the total extension over the gage 
length and the mean strain is computed by dividing the extension 
by the gage length. If the variation of strain along the gage length 
is not nearly linear, the measured value may differ considerably from 
the strain at the middle of the gage length. This difference may be 
reduced by the use of a shorter gage length but as a result the sensi- 
tivity of the gage would be decreased and the gage would be more 
difficult to construct. 

Another example is the measurement of heat capacity; the heat 
capacity, dQjdt, of a body at temperature t is usually approximated 
by AQ/Ai, where AQ is the quantity of heat required to change the 
temperature of the body by an amount At, where the interval At 
includes t. 

A method of correcting such data w r as obtained by Strutt [l]. 1 
Strutt's formula is equivalent to the first correction-term of eq 
12 in this paper. It will be seen that this is inadequate in many 



1 Figures in brackets indicate the literature references at the end of this paper. 
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cases. An analogous method for cases where the function is defined 
for discrete values of the argument only is King's formula for quin- 
quennial sums [2]. 

Runge [3] has proposed a method of correcting observed spectral- 
energy distributions which has been applied by Stang [4], and by 
Runge to some data of Paschen. In this case the observed intensity 
is a double integral of the desired intensity, the error arising from the 
finite width of both the collimating slit and the bolometer strip. 
Runge's formula is not applicable to the problem considered here. 

The statistical data which specify frequency distributions are often 
obtained in an approximate form analogous to the cases of strain and 
heat capacity. Thus, in a table of heights of individuals, the entry 
corresponding to 70 in. might be the number of individuals whose 
heights lie between 69K and 70K in. A method of computing the 
"true" moments of the distribution from the "rough" moments ob- 
tained from the data has been given by Sheppard [5]. Sheppard's 
work is the closest approach to a solution of the present problem found 
in the literature, but the method is in general not directly applicable. 

II. DERIVATION OF THE FORMULAS 

1. CENTRAL-DIFFERENCE FORMULA 

Suppose that it is desired to determine the values of f(x) correspond- 
ing to successive values of x differing by w, and that the method of 
measurement is such that the result, F{x), is the average value of f(x) 
over the interval g from x—g/2 to x-\-g]2. Then 

kw 

f^LVJ J {x)dx ' (i ) 

JX — ; 



kw 

' 2 



where k=g/uK This may be written symbolically as 

JW-is*^ (2) 

where D~d/dx and 5 k is defined from 

w««)-<(*4t)-*( , -t} (3) 

Eq 1 may be written 

j{x)=kwyF{x). (4) 

It is desired to calculate/Or) in terms of F(x) and its successive central 
differences, d n F(x), defined by repeated application of 

8 n iP(x)=d[d n - 1 ^(x)] 

to eq 3 with fc=l. The operators 5 k and D in eq 4 are thus to be re- 
placed by equivalent expressions in terms of the operator 6. These 
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expressions are obtained by expanding the right-hand side of eq 3 in 
T ardor's series about x, thus: 

kw kw 

or 

5*==2sinh^Z), (5) 

from which, with fc=l, 

D=- sinh" 1 J (6) 

w 2 w 

The substitution of Z) from eq 6 into eq 5 gives 

5*=2 sdnh (k mrihr 1 1\ (7) 

Replacement of D and 5* in eq 4 by the values given in eq 6 and 7 
gives 

k sinh -1 - 

/(*)=- -7 --s\ F W- (8) 

sinh( k sinh -1 -^] 

This is in the form f(x) =p csch^> F(x), where p=k sinh _1 (5/2), and by 
expansion of p csch p in MaclaurhVs series may be written 

f(x) =[l +2± y i) *®*~*~ 1) B in - 1 (k smh-'-Q 2 "]^), (9) 

where B 2n -i are the successive Bernoulli's numbers. 2 Expansion of 
sinh _1 (5/2) in Maclaurin's series in eq 9 gives the double series 

[l( 1 +g ( - 1)m (2^Pgn + D (5/2)2m )] J "K ) ' (10) 

or 

fM-n Py. gpO+W^ P(448+196**+31fc«) "1 

j W [_ A 24° + 5760 d 967,680 d +> - J W ' (Ll) 

When fc=l, which is the case when the intervals g are contiguous, 
eq 11 reduces to 

Kz) =(l- 2 V 2 +6i0 54 -7ll8 60 + • • "W (12) 

The numerical values of the coefficients of the first three correction 
terms in eq 11 are given for four values of k, in table 1. 

2 What are hero denoted by B 2n -i are sometimes denoted by B n . 
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Table 1.— Values of C 2 , C iy and C 6 infix) = (1 + C 2 5 2 +C 4 5 4 + C fi 5 6 +. . . ) F(x) 



k 


-c 2 


Q 


-c 6 


k 


-c 2 


d 


-c« 


1 


0.0417 
.167 


0. 00469 
.0333 


0. 000714 
. 00714 


3 

4 


0. 375 
.667 


0.130 
.367 


0.0439 


2 


.191 







2.1DESCENDING-DIFFERENCE FORMULA 

Equations 11 and 12 cannot be used to calculate j(x) at x=b, 
where b is near a boundary of the interval of x for which/Oc) is defined, 
or where values of x on only one side of b are accessible to the method of 
measurement, because the necessary central differences are not avail- 
able. For such cases an ascending- or descending-difference formula 
may be used. It is convenient to consider the measured value F(x) 
to correspond to the initial point of the interval g=kw, so that 



or, symbolically, 



1 fz+kw 



JPM— 1 A /(*) 



(13) 



(14) 



where D = d/dx and A*, the descending difference, is defined from 

A k f(x) = i(x+kw)-4<(x). (15) 

Equation 14 may be written 



f(x)=kw^-F(x)- 



(16) 



It is desired to calculate j{x) in terms of F(x) and its successive de- 
scending differences A n F(x), defined by repeated application of 

A n $(x) = A[A n - 1 f(x)] 

to eq 15 with 4=1. The operators A k and I) in eq 16 are thus to be 
replaced by equivalent expressions in terms of the operator A. Equa- 
tion 15 may be written 



or 



and 



(l + A k )t(x) = t(x+kw)=e k » D t(x) f 
l+A k =e kwD , 



1+A=e w2> . 

Equation 17 and 18 give 

wD=log (1 + A), 



and 



A t =(l + A)*-1. 



(17) 
(18) 

(19) 
(20) 
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Replacement of D and A k in eq 16 by the values given in eq 19 and 
20 gives 



JW (1-1- A)* yF\ X ) m 

V 



(21) 



This is in the form/(^)= p _. F{x), where p=k log (1+A), and by 
expansion of f]{eP — 1) in Maclaurin's series may be written 

/(x)=jl-|log (l + A)+g^g^[*log (i +A)P JF( x ), (22) 

where B 2n -\ are the successive Bernoulli's numbers. 3 Expansion of 
log (1 + A) in Taylor's series in eq 22 gives the double series 



(2n)! 



Hs(-i) m+1 ^ \F(x), 



or 



(23) 



^)=[ 1 -2 A +|( 1 +|> 2 -6( 1 +l> 3 +|( 1 +lV"-|> 4 
_k/l 5jc_J^\ Jc/ lJ7k_UP J^\ 1 

2\.5 + 36 180 )* + 12V f 180 360 +2520^ + " " J W ' ( ^ 4) 

When 4=1, which is the case when the intervals g are contiguous, 
eq 24 reduces to 



y W -(l-jA+|A»-iA"+jA*-jA»+jA- 



.)F(x). (25) 



The numerical values of the coefficients of the first six correction 
terms in eq 24 are given, for four values of k, in table 2. 



Table 2 


.< — Values of Ci, 


c 2 , 


etc. in 


/(aO = (l+C,A + C 2 A 3 +C 8 A 8 + • 


. .)F(x) 


k 


-Ci 


d 


-d 


d 


-C« 


a 


1 






0.500 
1.000 
1.500 
2. 000 


0.333 

.834 

1.500 

2.333 


0.250 

.667 

1.250 

2.000 


0.200 
.533 
.950 

1.367 


0.167 
.433 
.700 
.800 


0.143 


2 


.360 


3... 


.526 


4 . 


.476 







III. ACCURACY AND LIMITATIONS 

The derivations of eq 11 and 24 involve the expansion of operators 
in terms of ascending powers (orders) of the symbols 5 and A, respec- 
tively. This procedure is strictly valid if the function Fix) upon 
which the expanded expression operates is a polynomial, because 8 
and A obey the commutative, associative, and distributive laws; and 
because, since the operation 5 or A on a polynomial of finite degree n 
reduces its degree by one, all differences of Fix) of order higher than 
the nt\\ vanish. The series therefore terminates, and is exact. 

However, if Fix) is not a polynomial, the formulas 11 and 24 can be 
applied only on the assumption that Fix) may be sufficiently well 
represented by a polynomial of degree n. The error involved in such 
a procedure might be estimated by considering the remainder after 

3 See footnote 2. 
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the nth. difference term of the series, but this is not possible when, as 
is usually the case, the characteristics of Fix) are unknown. The 
justification for the use of formulas 11 and 24 therefore rests on the 
fact that distributions encountered in practice are frequently such 
that they may be approximated to by polynomials of low degree. 

It is obvious that, so far as the data F(:r) are concerned, the function 
f(x) is undetermined to the extent of an additive arbitrary periodic 
function of period g the average value of which over any interval g 
is zero. This indetermination can be reduced only by reducing the 
interval g over which /(#) is averaged in the measurement. Of course, 
a lower limit to g is set by other considerations. 

The value for the interval w between observations may ordinarily 
be the same as that which would for various reasons have been used 
if the corrections were not to be applied. In cases where there are 
not enough data to provide differences of sufficiently high order, w 
can be decreased. The highest order of differences which can be 
advantageously used is determined by the accuracy of the data. The 
successive differences are increasingly affected by the errors of meas- 
urement, and the differences of a certain order, and all succeeding 
differences, will consist mostly of accumulated error. It may be 
desirable to graduate the data F{x) before applying eq 11 or 24. 

IV. EXAMPLES 

It was thought desirable to include the following two examples 
which indicate roughly the order of accuracy to be expected from the 
use of eq 11 and 24. For these examples hypothetical data were 
obtained by computing mean values of strain over intervals of length 
g from theoretically known strain distributions. These mean values 
represent data that would be obtained by a perfect strain gage of 
gage length g. 

1. CENTRAL-DIFFERENCE FORMULA 

Consider (fig. 1) a very long isotropic elastic strip of breadth 2b 
containing a centrally located circular hole of diameter b. The strip 
is in a state of generalized plane stress such that the stress at cross 
sections remote from the hole is of magnitude xx —>T f uniform across 
the section and normal to it. The boundaries of the strip are free. 
The curve of figure 1 shows according to Howland [6] the variation of 
stress along an edge of the strip near the hole. Suppose that in a 
certain specimen the stresses along an edge were actually as given by 
the curve of figure 1 and it were desired to determine them by means 
of a strain gage. (Since the edge of the strip is a free boundary, the 
stress is proportional to the strain.) The data that would be obtained 
by the use of a strain gage of gage length (2/3)6 and the results of 
applying eq 11 are shown in figure 1. It is evident that the corrected 
values conform much more closely than the data to the original 
distribution. 

2. DESCENDING-DIFFERENCE FORMULA 

Consider (fig. 2) an isotropic elastic flat circular plate of radius a 
and thickness h<£.a, clamped around the edge and subjected to a 
hydrostatic pressure p = 18. %E(}i A /a 4 ) on one face. (E is Young's 
modulus of elasticity and Poisson's ratio is taken as 1/4.) The curve 
of figure 2 shows according to Nadai [7] the variation along a radius of 
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.given by strain gage, gage length 2/3 b 



4- 



4- 



d corrected by eq II with A = 2 

I I 

* corrected by eq II with k = 4 




? --- 
5 h 



<i 



h 2 



12a/ 



Figure 1. — Problem of a tension member containing a hole. 

i where e r is the normal strain in the unloaded face in the radial 

direction. Suppose that in a certain specimen the radial strains were 
actually as given by the curve of figure 2 and it were desired to deter- 
mine them by means of a strain gage. If the value of the strain at the 
clamped edge were desired, eq 24 would be used with the trans- 
formation x=\—u. The data that would be obtained by the use 
of strain gages of gage lengths (2/5)a and (4/5)<x and the results of 
applying eq 24 are shown in figure 2. 

The corrected values conform much more closely than the data to 

the original distribution. Furthermore, the value ofe r \j7y-2 at the 

clamped edge, which is inaccessible to the strain gage, is accurately 
obtained. 

3. NOTATION 

Various notations for differences are in use. Tables 3 and 4 show 
how the differences used in this paper are computed. The tables are 
arranged so that the difference corresponding to any entry, F(x), is 
on a horizontal line through that entry. 

The interval between successive values of the argument x is constant. 
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-100 



Figure 2. — Problem of a circular flat plate subjected to hydrostatic pressure. 
Table 3. — Central differences 



Fix) 


8Fix) 


&F(x) 


&F(x) 


8*Fix) 


F(xo) 


6F(xi,2) =Fix l )-Fixo) 








F(xO 


dFixw) <=Fix 2 )-Fixi) 


82FixO=8Fix m ) -5Fix l/2 ) 


SSFixw) =8*F(Z2)S*F(xi) 




F(x 2 ) 


SFix b / 2 ) =Fixz)-Fix 2 ) 


82Fix 2 )=8Fixw) -6F(zm) 


8*Fixm) =8 i Fix 3 )-82Fix2) 


8 i Fix 2 )^Fix 5 /2) -&F(zz/i) 


F(xz) 


8Fix m ) =FiXi)-FiXi) 


8*Fix z )=8Fixv2) -8Fix b / 2 ) 


8*F(xy2) =5 2 FiXiy82Fix3) 


8*Fixa)=83Fix 7 /2) -8^Fix 5 /i) 


F{x,) 


BFix m ) =Fix b )-FiXi) 


d*Ffa)^F(xw) -8Fix V 2) 


8*F(xw) =8 2 Fixt,)-8iFiXi) 


VFtefi'&Fixw) -8 3 Fix 7/i ) 


F{x,) 


8F(zu/2)=F(xs)-F(zi) 


82Fix 5 )=8Fixn/2)-8Fix V2 ) 


8*Fixn/2)=8 2 Fix6)-82Fix!i) 


8*F(X5)=8SFixn/2)-8SFix 9 /2) 


Fix,) 


SFixiz/2)=Fix 7 )-Fix 6 ) 


82Fix*)=SFix W 2)-8Fix W 2) 






Fixi) 
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Fix) 


AFix) 


A*Fix) 


A*Fix) 


A*Fix) 


F(x») 


AFixo)=Fix l )-Fix Q ) 


A*Fixo)=AFix{)-AFiXo) 


AtFix )=A*Fixi)-AiFix ) 


A*Fixo)=A2Fixi)-ASFiX()) 


FCxi) 


AFixi)=Fix 2 )-Fixi) 


AtFixi)=AFix 2 )-AFixi) 


A^Fix^A^FiXihA^Fixi) 


A*Fixi)=A3Fix 2 )-ASFixi) 


Flu) 


AFix 2 )=Fix3)-Fix 2 ) 


AiFix 2 )=AFixi)-AFiX2) 


A3F(x 2 )=A2Fix-i)-A , >FiX2) 


A*Fix 2 )=AU<Xx3)-A* F(Xi ) 


F(x 3 ) 


AFix 3 )=FiXi)-FiX3) 


AiFixz)=AFiXi)-AFiXi) 


AZFix^AtFixJ-AtFiX's) 


A^Fix 3 )=A^Fix i )-AU<\.i: i ) 


F(Xi) 


AFix^FixzhFiXt) 


A2FiXi)=AFix 5 )-AFiXi) 


AAFix^AZFixahAtFix^ 




F(x,) 


AFix b )=Fix b )-Fix b ) 


AiFiX!,)=AFix 6 )-AFix 5 ) 






Fixt) 


AFix^Fixd-Fixt) 








Fix 7 ) 
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